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We report results on the ratio of mid-rapidity anti-proton to proton yields in Au+Au collisions 
at ^/snn = 130 GeV per nucleon pair as measured by the STAR experiment at RHIC. Within 
the rapidity and transverse momentum range of \y\ < 0.5 and 0.4 < p t < 1.0 GeV/c, the ratio 
is essentially independent of either transverse momentum or rapidity, with an average of 0.65 ± 
0-01( s t a t.) ±0.07( syBt .) for minimum bias collisions. Within errors, no strong centrality dependence is 
observed. The results indicate that at this RHIC energy, although the p-p pair production becomes 
important at mid-rapidity, a significant excess of baryons over anti-baryons is still present. 

PACS numbers: 25.75.Ld 



Lattice Quantum Chromodynamics calculations pre- 
dict that at sufficiently high energy density a phase tran- 
sition from hadronic matter to a state of deconfincd 
quarks and gluons fj], |) will occur. To create and study 
this deconfined state is the primary goal of the heavy-ion 
collision program at the Relativistic Heavy Ion Collider 
(RHIC) Q. However, the formation of such a deconfincd 
state depends on the initial conditions of the matter cre- 
ated at the early stage of heavy- ion collisions. Baryon 
number transport (or stopping), achieved mostly at the 
early stage, is one of the important observables j^, [I| for 
high-energy collisions, since the degree of baryon stop- 
ping affects the overall dynamical evolution of these col- 
lisions. It affects the processes of initial parton equilibra- 
tion ||, 0, ||, particle production thermal and/or 
chemical equilibration j| , and the development of collec- 
tive expansion jl0|, [ll] . Information on baryon transport 
may be experimentally accessed by the measurement of 
the ratio of the anti-proton to proton yields (p/p)- 

In this letter, we report results on the inclusive p/p 



ratio in Au+Au collisions at the center of mass energy 
y / snn = 130 GeV per nucleon pair measured by the 
Solenoidal Tracker at RHIC (STAR). 

The STAR detector consists of several detector sub- 
systems in a large solenoidal magnet. The STAR Time 
Projection Chamber (TPC) Jl2| is 4.2 meters long, with 
a 50 cm inner radius and a 2 meter outer radius. For 
minimum ionizing particles in P-10 (90% argon, 10% 
methane) gas, approximately 45 primary electrons are 
produced per centimeter of track. This ionization was 
measured on 45 pad rows. The system could accommo- 
date ionization from 200 MeV/c protons without satura- 
tion and its noise level was about 5% of minimum ion- 
izing. More details of the detector can be found else- 
where jl3[ |lj]. For the data taken in the year 2000 run 
and presented here, the main setup consists of the TPC, 
a scintillator trigger barrel (CTB) surrounding the TPC, 
and two Zero Degree Calorimeters |l5| located up- and 
down-stream along the axis of the TPC and beams. The 
TPC was operated in a 0.25 Tesla magnetic field. The 
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CTB measures the total energy deposition in the scin- 
tillator from charged particles, and the ZDCs measure 
beam- like neutrons from fragmentation and/or evapora- 
tion of the colliding nuclei. The coincidence of the ZDCs 
and RHIC beam crossing was the experimental trigger 
for the minimum bias events used in this analysis. 

The collision centrality was determined off-line from 
the measured total charged particle multiplicity in the 
pseudorapidity range of \rj\ < 0.75. As was done in |pl| , 
the total charged multiplicity was scaled by the maxi- 
mum value measured, and the distribution of the scaled 
multiplicity was subdivided into eight centrality bins. 

For this analysis, 68k minimum bias events were used 
with an event vertex \z\ < 50 cm. Protons and anti- 
protons were selected according to specific energy loss 
(dE/dx) in the TPC up to a momentum of 1 GeV/c. 
A mean dE/dx for each track was obtained by aver- 
aging the lower 70% of the measured dE/dx values. 
This selection reduces our sensitivity to large fluctua- 
tions in the dE/dx measurements. At a momentum of 
0.5 GeV/c, the dE/dx resolution was found to be 8% 
for a typical long track in the TPC. Figure [j] shows 
the mid-rapidity (\y\ < 0.1) particle raw yields as a 
function of Z = log[(dE/dx)Exp/{dE /<1x)bb\, where 
(dE/dx)sB is the expected dE/dx value. The proton 
mass was used in the Bethe-Bloch (dE/dx)BB calcula- 
tions so the Z-distributions of proton (anti-proton) are 
centered about zero in those plots. The dashed- lines are 
the multi-Gaussian fits to the Z-distributions (including 
pions, kaons, and protons). Protons and anti-protons 
with p t > 175 MeV/c were reconstructed. Within the ac- 
ceptance the typical momentum resolution was 2% and 
the systematic uncertainty in the proton (anti-proton) 
yields, extracted from the dE/dx fitting, was less than 
2%. Tracks were selected if their distance of closest ap- 
proach to the primary vertex were less than 3 cm, and if 
they had at least 15 out of 45 possible TPC space-points. 
Several tests showed that the p/p results were not sensi- 
tive to variations of these cuts within reasonable ranges. 
The raw yields of protons and anti-protons were obtained 
by fitting the Z distributions (see Figure 1) for each given 
(y,Pt) bin, and were used to obtain the p/p ratio. 

The STAR TPC is symmetric about mid-rapidity and 
has full azimuthal coverage. Due to this symmetry, most 
of the detector effects are the same for protons and anti- 
protons and cancel in the p/p ratio. However, there are 
two important differences which are mostly dominant at 
low momentum: (i) background protons are produced in 
the detector materials via hadronic interactions [fl6[ , and 
(ii) some of the anti-protons are absorbed in the detector 
materials. 

Background protons are evident from the distribution 
of the distance of closest approach (dca). The distribu- 
tion is peaked at small dca and has a flat tail from sec- 
ondary production which extends in to the peak region 
where primary tracks are centered. Since the background 
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FIG. 1: Mid-rapidity particle yields as a function of Z — 
log[(dE /dx)Exp/(dE /dx)BB\ from minimum bias Au+Au 
events. Top two plots (a)-(b) and bottom plots (c)-(d) are for 
positive particles and negative particles, respectively. Pro- 
ton mass was used in the Bethe-Bloch (dE / dx)BB calcula- 
tions so the Z-distributions of proton are centered about zero. 
Dashed-lines are the multi-Gaussian fits to the measured dis- 
tributions. 

shape of the dca distribution at small dca (< 3 cm) is not 
accessible from the data, a Monte Carlo (MC) simulation 
(GEANT |0) was employed using as input particles gen- 
erated by both HIJING |§] and RQMD || models. An 
empirical form was found for the background: 



b ackground (dca ) 



1 + exp 



dca — dcao 



(1) 



where the parameters dcag and a were obtained by fit- 
ting to the MC results. In general these parameters are 
momentum dependent and the normalization factor was 
obtained from the data. In the momentum region around 
0.2 GeV/c, the number of background protons was found 
to be twice the number of primary protons. This leads to 
a rather large systematic uncertainty in the raw proton 
yield. Therefore, in this analysis we limited ourselves to 
the region where systematic errors are below 10% lead- 
ing to a lower p t cut-off at 0.4 GeV/c. At the high mo- 
mentum end, p w 1 GeV/c, the dE/dx method becomes 
insufficient for particle identification. 

The anti-proton absorption loss was estimated via: 



absorption = 1 - cxp(-(T anni/ 9 t p/p t ), 



(2) 



where pt is the transverse area density of nucleons in 
the materials, p and p t are the anti-proton total mo- 
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mentum and transverse momentum, respectively. The 
annihilation cross section was parameterized as cr a nni = 
1.2 CTtot/v^j where ^/s is the center-of-mass energy of 
the p-p pair in GeV, and a power-law parameterization 
for the total p-p cross-section <7 t ot from Ref. [^0| was used. 
For consistency, we also checked the above parameteriza- 
tion for anti-proton absorption loss in the TPC with the 
results of a full MC simulation for the STAR TPC. The 
results are consistent and within the kinematic region of 
0.4 < pt < 1.0 GeV/c, the p absorption loss was found 
to be less than 5%. 

In this paper, we report the inclusive p/p ratios and no 
attempts have been made to correct for feed-down pro- 
tons and anti-protons from hyperon weak decays. Such 
corrections would inevitably depend on the assumptions 
for hyperon and anti-hyperon yields and momentum dis- 
tributions. However, if anti-hyperon to hyperon ratios 
are the same as the p/p ratio, then weak-decay feed- 
downs would not affect the p/p ratio. Using the HI- 
JING 1 18 results as input, MC simulations show that 
about 81% of the A-decay protons were reconstructed 
with dca < 3 cm and 65% with dca < 1 cm. As ex- 
pected, the fractions are the same for anti-protons. For a 
systematic check, some of the hyperon feed-down effects 
were studied from the data by varying the dca cut. The 
results show that the p/p ratio decreases by about 2% 
from dca = 3 cm to 1 cm cut. These effects were not 
included in the systematic errors estimated below. 

For this measurement, systematic errors mainly come 
from two sources: (i) The systematic uncertainty in pro- 
ton background in relatively low transverse momentum 
region (p t < 0.4 GeV/c). As mentioned above, the back- 
ground protons were estimated via the detailed GEANT 
jlTj simulation of the TPC. By varying the dca cut, the 
systematic error on the determination of number of pro- 
tons is < 10% at p t ~ 0.4 GeV/c and drops to < 2% 
at p t ~ 0.6 GeV/c. (ii) As one will see in Figure 2(b), 
there is an asymmetry between positive and negative ra- 
pidities. The asymmetry is less than 7%. The origin of 
the asymmetry is not fully understood and therefore it is 
included in the overall systematic errors for the p/p ratio. 
Other uncertainties like contamination in particle identi- 
fication (most at high p t ) and the centrality dependence 
of the systematic errors are estimated much smaller, and 
are not included in the overall systematic errors. 

Figure ||(a) shows the minimum bias p/p ratio as a 
function of p t in the rapidity range of 0.3 < \y\ < 0.4. 
The ratio is 0.65 and is consistent with a constant value 
in the pt range of 0.4 < p t < 1.0 GeV/c. The systematic 
error is estimated to be 10% for the lowest p t bin, mainly 
due to proton background. At higher p t the systematic 
error is less than 7%. Figure ||(b) shows the p/p ratio 
integrated over 0.6 < pt < 0.8 GeV/c as a function of 
rapidity. In this kinematic region, systematic errors are 
estimated to be less than 10%. Within errors, the ratio 
is consistent with a constant in the measured rapidity 



range. 

Figure |^(c) shows the centrality dependence of p/p in- 
tegrated over \y\ < 0.3 and 0.6 < p t < 0.8 GeV/c. Al- 
though the ratio is consistent with a constant, there is an 
indication of a systematic drop from peripheral to cen- 
tral collisions. At lower bombarding energies, the value 
of the p/p ratio decreases by a factor of 2 with increas- 
ing centrality for Au+Au collisions at the AGS [pi], |22] ] 
and a factor of 1.6 for Pb+Pb collisions at the SPS 
p3| , |2~i| , |25| , p6| ^7j. These low energy results are con- 
sistent with more baryon stopping and/or nucleon-p an- 
nihilation in central collisions relative to peripheral col- 
lisions. A similar picture may also apply to the RHIC 
results presented here. 
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FIG. 2: The anti-proton to proton ratios (a) as a function 
of transverse momentum pt over the rapidity range 0.3 < 
\y\ < 0.4; (b) as a function of rapidity within 0.6 < pt < 
0.8 GeV/c; and (c) as a function of the collision centrality 
within (\y\ < 0.3 and 0.6 < p t < 0.8 GeV/c). (a) and (b) are 
for minimum bias collisions. Errors are statistical only. The 
overall systematic errors are estimated to be 10%. 

At this RHIC energy, the p/p ratio is significantly 
smaller than unity over the measured centrality range, 
indicating an overall excess of protons over anti-protons 
in the mid-rapidity region. This implies that a certain 
fraction of the baryon number is transported from the 
incoming nucleus at beam rapidity to the mid-rapidity re- 
gion even in peripheral Au+Au collisions at ^/snn = 130 
GeV. Thus at this RHIC energy, the mid-rapidity region 
is not yet net-baryon free. 

On the other hand, there is a dramatic increase in the 
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FIG. 3: Mid-rapidity anti-proton to proton ratio (p/p) mea- 
sured in central heavy-ion collisions (filled symbols) and el- 
ementary p + p collisions (open symbols). The left end of 
the abscissa is the p-p pair production threshold in p + p 
(V S NN = 3.75 GeV). The RHIC data is the most central 
point from Fig. 2c. Errors, either shown or smaller than 
the symbol size, are statistical and systematic errors added in 
quadrature. 



mid-rapidity p/p ratio in central heavy ion collisions in 
going from AGS folp =0.00025±10%) |H) to SPS {p/p w 
0.07±10%) H, §7) and to RHIC. This is demonstrated 
in Fig. ^ where the central heavy-ion results are shown 
as a function of the center of mass energy ^/snn- For 
comparison, also shown in Fig. |], are the p/p ratios in 
p + p collisions at mid-rapidity and averaged over 0.35 < 
p t < 0.6 GeV/c ]2S|| . The kinematic coverage for the 
p + p collisions is similar to this measurement. The value 
of the p/p ratio from RHIC is close to or even larger than 
that in p+p collisions at v / snn = 53 GeV. Note that the 
beam rapidities at ^/snn = 53 and 130 GeV are 4.0 and 
4.9, respectively. 

For comparison, the HIJING(vl.35) model p8] pre- 
dicts a p/p ratio of approximately 0.8 for Au + Au cen- 
tral collisions at ^snn = 130 GeV. If the baryon junc- 
tion mechanism [Q is considered, the ratio is 0.75. In 
contrast, the Relativistic Quantum Molecular Dynam- 
ics RQMD(v2.4) model [|l9| predicts a p/p ratio increas- 
ing with p t with an average value of approximately 0.5. 
Note that the HIJING model is motivated by perturba- 
tive QCD but does not have late stage rescattering, and 
the RQMD model has late stage hadronic rescattering. 
As a result, the RQMD calculations predict a pt depen- 
dence of the p/p ratio (from 0.2 to 0.55 in p t < 1 GeV/c) 
while a constant ratio is observed from the HIJING cal- 
culations. 

In summary, we have reported the ratio of the mid- 
rapidity anti-proton to proton yields in Au+Au collisions 
at y^SNN = 130 GeV measured by the STAR experiment. 
Within the rapidity and transverse momentum range of 



\y\ < 0.5 and 0.4 < p t < 1.0 GeV/c, the ratio is es- 
sentially independent of either transverse momentum or 
rapidity. In this kinematic range, the average p/p ratio is 
found to be 0.65±0.01( stat )±0.07( syst .) for minimum bias 
collisions. No strong centrality dependence is observed in 
the centrality range measured. The value of the p/p ratio 
indicates that although p-p pair production becomes im- 
portant at mid-rapidity, a significant excess of baryons 
over anti-baryons is still present in heavy ion collisions 
at RHIC. Comparisons of our result to heavy ion results 
at lower energies indicate that the mid-rapidity p/p ra- 
tio in heavy ion collisions increases significantly with the 
collision energy. 
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